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I.   Introduction 

This  report  describes  position  estimation  procedures  that 
are  based  on  models  which  relate  positional  uncertainty  to 
various  measurement  and  estimation  errors.   The  procedures  were 
developed  to  be  used  in  analyses  that  relate  the  effect  of  posi- 
tional uncertainty  on  tactical  performance  through  such  factors 
as  weapon  accuracy  as  well  as  to  be  used  operationally. 

In  the  models,  positions  are  on  a  plane  surface  (flat  earth 
model).   Because  of  this  condition,  the  models  are  not  intended 
for  use  in  situations  where  the  earth's  figure  is  significant. 
In  addition,  positional  errors  are  determined  by  independent 
normally  distributed  random  variables  with  known  means,  variances 
and  covariances.   The  support  for  this  condition,  other  than  its 
mathematical  convenience,  is  that  it  has  been  used  by  others, 
for  example,  see  Reference  1.   To  use  the  models,  one  is  required 
to  specify  the  means,  variances  and  covariances  of  the  error 
random  variables.   For  these  models,  this  can  be  done  by  speci- 
fying the  systematic  errors  (biases)  and  the  error  magnitudes 
(standard  deviations) . 

The  procedure  that  is  described  in  Section  II  relates  position 
estimates  that  are  based  on  bearings  on  or  from  stations  to 
station  position  uncertainty.   One  application  for  the  model  is 
the  analysis  of  the  effect  of  sonobuoy  position  uncertainty  on 
position  estimates  determined  with  passive  directional  sonobuoys . 

The  procedure  that  is  described  in  Section  III  relates  posi- 
tion estimates  that  are  based  on  lines  of  position  to  line  of 


position  uncertainty.   In  addition  to  estimating  the  uncertainty 
in  position  estimates  based  on  celestial  observations,  the  pro- 
cedure could  be  used  to  determine  error  ellipses  for  LORAN  fixes 
if  the  standard  deviation  values  required  by  the  procedure  could 
be  obtained. 

The  procedure  that  is  described  in  Section  IV  can  be  used 
to  combine  position  estimates  from  various  sources.   The  proce- 
dure which  is  based  on  conditions  that  should  not  be  too  restric- 
tive  in  most  cases  provides  both  a  composite  position  estimate 
and  error  ellipse. 


II.   Station  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this   section  relates 
position  estimates  that  are  based  on  bearings  on  or  from  stations 
to  station   position  uncertainty.   The  procedure  is  based  on  a 
model  that  is  an  extension  of  one  that  is  described  on  Reference 
2.   The  model  is  defined  as  follows:   Each  station  position  error 
is  determined  by  an  independent  bivariate  normal  distribution 
with  a  zero  mean  vector  and  a  known  covariance  matrix.   Observed 
bearing  lines  on  or  from  a  station  are  parallel  to  true  bearing 
lines.   The  distance  of  each  observed  bearing  line  from  its  true 
bearing  line  is  determined  by  an  independent  normal  distribution 
with  a  zero  mean  and  a  standard  deviation  a. 

Because  a  station's  position  error  is  determined  by  a  bivar- 
iate normal  distribution,  the  perpendicular  distance  s  between  a 
line  at  the  assumed  location  of  an  observed  bearing  line  and  the 
observed  bearing  line  is  determined  by  a  normal  distribution. 
The  relation  between  the  bivariate  normal  distribution  that 
describes  the  station's  position  uncertainty  is  indicated  in 
Figure  1.   In  the  figure,  the  positive  y-axis  direction  is  north, 
the  positive  x-axis  direction  is  east,  and  the  origin  of  the 
coordinate  system  is  at  the  assumed  station  position.   The  x1-'- 
coordinate  system  is  oriented  so  that  the  positive  x'axis  is 
coincident  with  the  major  axes  of  the  elliptical  contours  of 
the  bivariate  normal  distribution  that  determines  the  station 
position  error  and  so  that  the  bearing  6  of  the  positive  x'-axis 
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Figure  1.   The  geometry  associated  with  the  determination  of  os, 
the  standard  deviation  of  the  normal  distribution 
that  determines  the  distance  between  the  assumed 
location  of  an  observed  bearing  line  and  the  observed 
bearing  line.   The  assumed  (mean)  position  of  the 
station  is  at  the  origin.   The  ellipse  represents  a 
contour  on  the  probability  density  surface  of  the 
bivariate  normal  distribution  that  determines  the 
station's  position. 


satisfies  the  condition:  0°    <$  <  180°.   The  anqle  0  is  an 
observed  bearing.   The  x "y "-coordinate  system  is  oriented  so 
that  the  positive  y"-axis  is  in  the  direction  of  the  observed 
bearing  line.   As  a  consequence  of  these  relationships, 

2       2     2  2     2 

a   =  G  .   sin  (6-9)  +  a.   cos  (6-q) 

S        X  y 

2  2 

where  a  .   and  a  .   are  the  elements  of  the  station  position 
x '        y 

error  covariance  matrix  relative  to  the  x ' y ' -coordinate  system. 

A  procedure  for  determining  position  estimates  that  are 
based  on  bearings  on  or  from  stations  at  known  positions  is 
described  in  Appendix  1  of  Reference  3.   The  procedure  is  based 
on  a  model  that  is  equivalent  to  one  that  is  described  in  Ref- 
erence 2.   The  procedure  that  is  described  in  this  section  is 
based  on  a  model  that  is  an  extension  of  it. 

The  procedure  in  Appendix  1  of  Reference  3  is  based  on  a 
model  in  which  a  station's  bearing  error  is  determined  by  a 
normal  distribution  with  zero  mean  (bias)  and  standard  devia- 
tion e.   The  bearing  error  is  related  to  the  distance  on  a  cir- 
cular arc  between  a  station's  true  bearing  line  and  observed 
bearing  line.   The  arc  is  on  the  circle  with  its  center  at  the 
station  that  passes  through  an  initial  estimate  of  an  object's 
position.   This  distance  is  determined  by  a  normal  distribution 
with  mean  zero  and  standard  deviation  a  =  re  where  r  is  the 
range  of  the  initial  estimate  from  the  station  and  the  standard 
deviation  e  is  measured  in  radians.   In  the  model,  arc  distance 


is  approximated  using  a  first  order  approximation  which  in 
effect  replaces  the  circle  with  its  tangent  line  at  the  initial 
estimate.   As  a  consequence,  the  distance  u  on  the  tangent  line 
between  the  observed  bearing  line  and  the  true  bearing  line  is 
determined  by  a  normal  random  variable  with  mean  zero  and  stan- 
dard deviation  a.   The  distance  u  can  be  expressed  in  terms  of 
w,  the  distance  on  the  tangent  line  between  the  observed  bearing 
line  and  the  initial  estimate  that  is  also  determined  by  a  normal 
random  variable  with  standard  deviation  a,  and  v,  the  distance 
on  the  tangent  line  between  the  true  bearing  line  and  the  initial 
estimate.   And,  as  shown  in  Appendix  1  of  Reference  3,  this 
distance  can  be  expressed  in  terms  of  the  unknown  coordinates 
of  the  object's  position. 

The  effect  of  station  position  uncertainty  is  accounted  for 
by  the  distance  s  between  the  observed  bearing  line  and  the 
assumed  observed  bearing  line.   In  the  model  s  is  determined  by 

a  normal  distribution  with  mean  zero  and  standard  deviation  a 

s 

as  given  above.   This  approximation  is  consistent  with  the  first 
order  approximation  of  arc  distance.   As  a  consequence  of  these 
two  approximations,  all  of  the  bearing  lines  are  replaced  by 
lines  parallel  to  the  line  joining  the  initial  estimate's 
position  and  the  station's  assumed  position,  both  of  which  are 
known  positions.   The  geometry  involved  in  shown  in  Figure  2. 
The  modified  relationships  resulting  from  the  introduction  of 
station  position  uncertainty  are  shown  in  Figure  3.   The 
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Figure  2.   The  geometry  of  the  position  estimation  model.   The 
bearing  3  and  range  r  of  the  initial  estimate  from 
the  assumed  station  position  have  the  role  of  3 
and  r  in  Appendix  1  of  Reference  3. 
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Figure  3.   The  quantities  u,  v  and  w  correspond  to  the  quantities 
u,  v  and  w  in  Appendix  1  of  Reference  3.   The  auxiliary- 
quantities  u1  and  w"  are  defined  in  the  figure.   The 


9        9 

replacement  of  o^  by  a-  + 


,2  in  the  procedure  in 


Appendix  1  of  Reference  3  is  justified  by  noting  that 
with  this  substitution  u'  and  w1  are  equivalent  to  u 
and  w. 


important  thing  to  note  is  that  u'  =  u  -  s  is  determined  by  a 

normal  distribution  with  mean  zero  and  standard  deviation 

2      2  h 

[o   +o   ]2  and  that  otherwise  U*  is  equivalent  to  U  with 

respect  to  the  procedure  in  Appendix  1  of  Reference  3.   As  a 

consequence  of  this,  the  procedure  can  be  extended  to  include 

2      2    h 

station  uncertainty  by  replacing  o  by  [a   +  o   ]2  where  ever 

it  is  used.   In  this  case,  o  =  re  where  r  is  the  range  of  an 
initial  estimate  of  an  object's  position  from  the  assumed  station 
position  and  e  is  the  bearing  error  (standard  deviation)  in 
radians  of  the  bearings  associated  with  the  station. 


III.   Line  of  Position  Uncertainty  and  Position  Estimates 

The  procedure  that  is  described  in  this  section  relates 
position  estimates  that  are  based  on  lines  of  position  to  line 
of  position  uncertainty.   The  procedure  is  based  on  a  model  that 
is  defined  as  follows:   Lines  of  position  are  straight  lines. 
Observed  lines  of  position  are  parallel  to  true  lines  of  posi- 
tion.  The  distance  of  an  observed  line  of  position  from  a  true 
line  of  position  is  determined  by  an  independent  normally  dis- 
tributed random  variable   with  known  mean  and  standard  deviation 
Lines  of  position  are  specified  in  terms  of  a  rectangular  coor- 
dinate system  with  the  origin  at  a  reference  point  as  shown  in 
Figure  4.   For  celestial  navigation,  an  appropriate  choice  for 
the  reference  point  would  be  the  assumed  position.   Since  bear- 
ing lines  are  lines  of  position,  this  model  differs  from  the 
model  that  is  described  in  Section  II  only  in  terminology. 
However,  operationally  the  use  of  the  model  that  is  described 
in  this  section  differs  in  the  way  the  standard  deviation  of  the 
distance  of  the  line  from  a  true  line  is  determined.   The  stan- 
dard deviation  associated  with  each  line  of  position  must  be 
specified.   If  this  can  be  done,  the  procedure  can  be  used. 
As  an  example,  suppose  that  the  values  are  a-,  for  the  first 
line  of  position  and  a~  for  the  second  where  a,  >  a^     and 
that  the  lines  of  position  intersect  at  a  90°  angle.   In  this 
example,  the  minimum  area  confidence  (probability)  region  is 
an  ellipse  that  is  centered  on  the  estimated  position.   And  the 
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position 


initial  estimate 
and  final  estimate 


1st  line  of  position 


2nd  line  of  position 


position 


Figure  4.   The  geometry  when  only  two  lines  of  position  are  used 
with  the  procedure  in  Appendix  1  of  Reference  3.   The  lines  of 
position  correspond  to  bearing  lines  for  observed  bearings  ■■->    and 
6o  with  respect  to  the  procedure.   The  procedure  determines   the 
size  and  the  orientation  of  elliptical  confidence  regions  of  con- 
fidence p,  the  location  of  the  initial  estimate  and  the  location 
of  the  final  estimate.   The  locations  of  the  estimates  correspond 
when  only  two  lines  are  used.   If  the  lines  were  obtained  from 
sextant  observations,  the  assumed  position  should  be  the  reference 
position   in  which  case  the  lines  of  position  would  be  determined 
by  the  azimuth  anges  a,  and  a„  and  the  distances  p,  and  p~. 


11 


estimated  position   is  at  the  intersection  of  the  bearing  lines. 

2 
For  a  confidence  (probability)  of  containment  of  1  -  exp(-k  /2), 

the  major  axis  of  the  ellipse  is  coincident  with  the  first  bear- 
ing line  and  it  is  of  length  2ka, ,  the  minor  axis  of  the  ellipse 

is  coincident  with  the  second  bearing  line  and  it  is  of  length 

2 
2ko „    and  the  area  of  the  ellipse  is  irk  a  o   . 
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IV.   A  Composite  Position  Estimate 

The  procedure  that  is  described  in  this  section  is  for 
combining  position  estimates  for  an  object  from  independent 
sources.   It  is  based  on  the  following  model:   The  rectanqular 
coordinates  of  each  position  estimate  are  determined  by  an  inde- 
pendent bivariate  normal  distribution  whose  covariance  matrix 
is  known  but  whose  mean  vector  is  not  known.   The  components  of 
the  mean  vector  for  each  of  the  distributions  are  the  unknown 
coordinates  x  and  y  of  the  object.   This  model  implies  that  the 
natural  logarithm  of  the  likelihood  function  for  a  set  of  n 
estimates  can  be  expressed  as  follows: 

n  ,   -1 

log  L  =  K  -  1/2  I   (x.  -  x)  £.  (x.  -  x) 

i  =  l   1        L 

where  K  is  a  constant,  x.  is  an  estimate  vector  with  components  x. 

and  y.  ,  x  is  the  common  mean  vector  with  components  x  and  y, 


the  unknown  coordinates  of  the  object,  and  L-     is  the  covariance 

The  maximum  likeli- 


2      2 
matrix  with  elements  a.~  ,  a.-   and  tf-~~ 


hood  estimates  x  and  y  of  the  unknown  coordinates  x  and  y  are 
the  solutions  of  the  two  simultaneous  linear  equations  deter- 
mined by 


(log  L] 


3x 


x=x 

y=y 


=  o 


and 


?  (log  L) 
9y 


x=x 

y=y 
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The  equations  can  be  written  as 

A  x  +  B  y  =  D 
B  x  +  C  y  =  E 

and  their  solution  as 

CD  -  BE  AE  -  BD 

x  = 2  y  = 2 

AC  -  B  AC  -  B 

where  A  =  la.,    B  =  Eb.,  C  =  Ec .  ,  D  =  E(a.x.  +  b.y.  ), 
1         l         l  11     r*  i  ' 

y  9 

E  =  E(b.x.  +  c.y.),  a.  =  a.A  /d . ,  b.  =  -  a.--/d.,  c.  =  a.-  /d . , 
11     ri  '   i     ly    l    l       lxy   i    i     ix  '     i  ' 

2    2       2 
d.  =  a.-  a.~  -  a.~~   and  all  of  the  sums  are  for  i  from  1  to  n. 
l     ix   ly     lxy 

Since  they  are  linear  combinations  of  the  estimates  x.  and 
y. ,  the  estimates  x  and  y  are  determined  by  a  bivariate  normal 
distribution.   Consequently,  all  that  is  required  to  determine 

this  distribution  is  its  mean  vector  with  components  y -  and  y- 

2     2 
and  its  covariance  matrix  with  elements  a~  ,  a-   and  a~~.   The 

x    y       xy 

mean  vector  is  determined  by 

y-  =E{(CD  -  BE) /(AC  -  B2)}  =  x 
and 

y-  =E{(AE  -  BD)/(AC  -  B2)}   =  y 

And,  the  covariance  matrix  is  determined  by 

a-2  =  E{(CD  -  BE)  -  E (CD  -  BE)}2/ (AC  -  B2)2 

9  9  2  2 

=  {C  (F+I  +  2H)  -  2CB(H+K+M)  +  B  (G+J+2K)  }/ (AC  -  B  ) 
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a-2  =  E{(AE  -  BD)  -  E(AE  -  BD)}2/(AC  -  B2)2 

=  {A2(G+J  +  2N)  -  2AB(H  +  K+M)  +  B2 (F+I  +  2L)  }/ (AC  -  B2)2 

a--  =  E{[  (CD-BE) -E  (CD-BE)  ]•  [  (AE-BD ) -E  (AE-BD)  ]  }/(AC  -  B2)2 
xy 

=  {  (AC+B  ) (H+K+M)-CB(F+I+2L)-BA(G+J+2N) }/(AC  -  B2)2 

where  F  =  ^a.-2  ,  G  =  ^a.-2,  H  =  ^b.a.^2  ,  I  =  ZbA^2, 

2    2  2  2 

J  =  Zc.  a.-  ,  K  =  Zb.c.a.   ,  L  =  Za.b.a.--.  M  =  Z(a.c.+b.  )a.- 
i   ly  l  l  ly  '        11  ixy'         i  i   i  ;  ixy 

and  N  =  Zb.c.a.--  where  all  of  the  sums  are  for  i  =  1  to  n . 

Ill  Xy 

By  using  arguments  given  in  Appendix  1  of  Reference  3,  one 
can  show  that  the  axes  of  the  elliptical  confidence  regions 
associated  with  x  and  y  are  coincident  with  an  x ' y ' -coordinate 
system  where  the  transformation  from  the  xy-coordinate  system 

to   this  system  is  the  coordinate  axes  rotation  through  the 

2      2 
angle  y  defined  by  tan  2y  =  2a--/  (a-   -  a-  )  .   For  a  confidence 

p,  from  Reference  3,  the  minimum  area  confidence  region  is  an 

2 

ellipse  with  semi-axes  k  o^,  and  k  a-   ,  and  area  irk   a-ia-, 

x  y ;  x   y 

where  k  =  [-2  ln(l-p)]'2  , 

2      2     2  .         2     2 

a-,   =  a-   cos  y  -  2a~   cos  y    sin  \  +  a~   sin  y, 
x'      x  xy  y 

2      2     2  2     2 

a-,   =  a-   sin  y    +  2a^"  cos  y  sin  y    +    a~   cos  y 
y*      x  xy  y 

and  the  center  of  the  ellipse  is  at  the  point  (x,y) .   In  this 

coordinate  system,  a-,-,  =  0  . 
1  x  y 
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The  above  equations  can  be  used  to  specify  a  composite  posi- 
tion estimate  in  terms  of  the  location,  orientation  and  size  of 
an  elliptical  confidence  region  which  is  generally  the  form  in 

which  position  estimates  of  the  kind  that  are  being  considered 

2      2 
here  are  specified.   But,  since  values  of  a^   /  a"   and  a~~ 

r  x      y        xy 

are  required  for  each  of  the  n  estimates  that  is  being  combined, 
a  way  is  needed  for  determining  these  values  given  the  orienta- 
tion and  size  of  an  elliptical  confidence  region.   A  procedure 
to  do  this  when  the  orientation  is  given  in  terms  of  the  direc- 
tion 6  of  the  major  axis  and  the  size  is  given  in  terms  of  the 
lengths  SMJ  and  SMI  of  the  semi-major  and  semi-minor  axes  and 
the  confidence  p  is  described  next. 

By  using  an  xy-coordinate  system  in  which  the  positive 
y-axis  direction  is  north  and  the  positive  x-axis  direction  is 
east  and  with  the  convention  0°    6  <  180°,  the  dependence  of 
the  value  of  the  rotation  angle  y    and  of  the  order  relation 

between  oA.  and  aA,  on  the  value  of  the  major  axis  direction  6 

x '       y  J 

is  indicated  by  the  following  table: 


0°  <  6  <   45°:     Y  =  6         and  o- ,  >  a- , 

y '     x 

45°  <  6  <  135°:     y    =  <5  -   90°  and  a-,  >  a~ , 

x '     y ' 

135°  <  6  <  180°:     y  =  6  -  180°  and  a-,  >    a~, 
-                '  y '     x 

With  an  order  relation  and  a  value  for  p,  values  for  a~ ,  and 

a~ ,  can  be  determined  with  values  for  SMJ  and  SMI.   With  values 

y 

2      2 

for  o~,  ,  a~,  and  y.    values  for  a~   ,  a~  and  a~-  can  be  found 
x'     y                     x     y       xy 

from  the  following  equations: 
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2       2     2         2     2 
o~   =  o~    cos   Y  +  o~    sin   y, 
x      x1  y '         '  ' 

2       2     2         2     2 

o~   =o~    sin   y  +  o~,   cos   y 
y     x'  y 

and 

2       2 

a~~  =  (a~,   -  a-   )  sin  y  cos  y 
xy     y      x1 

2 
which  can  be  obtained  by  inverting  the  equations  above  for  o- , 

2 

and  o~  . 

Y 

As  an  example,  suppose  the  data  in  the  following  table 
represent  three  independent  position  estimates: 


k 
2 
2 

J 

Here,  distances  are  in  nautical  miles  and  p  =  .86  in  each  case 
For  this  example  with  values  in  square  nautical  miles: 

a,-2  =  264.53,  a  -2  =  159.42  and  a _—  =  98.89 
lx  '   ly  lxy 

0o-2  =  321.35,  a  -2  =   51.15  and  a  --  =  -78 
2y  2y  2xy 

a  -2  =  249.20,  o  -   =  389.30  and  a.--  =  -173.4 
3x  '   3y  3xy 

These  values  give  the  following  composite  estimate: 

x  =  -2.46  nautical  miles  and  y  =  12.26  nautical  miles. 


X 

Y 

6 

SMJ 

SMI 

1st 

-3.7 

18.1 

59° 

36 

20 

2nd 

11.8 

8.4 

105° 

37 

11 

3rd 

0 

0 

146° 

45 

23 
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For  this  case,  o~  =  16.92  nautical  miles,  o~   =  6.21  nautical 

x  Y 

miles  and  o~^  =  -55.1  square  nautical  miles.   And,  for  k  =  2 : 
xy  ^ 

SMJ  =  33.85  nautical  miles,  SMI  =  12.42  nautical  miles  and 
6  =  102°.   The  composite  confidence  region  and  its  three  com- 
ponent confidence  regions  are  shown  in  Figure  5. 

As  a  second  example,  suppose  each  position  estimate  is  deter- 
mined by  a  circular  normal  distribution.   Then  a--  =  a.~  =  a. 
■*  ix     ly     i 

and  a.--  =  0  for  i  =  1  to  n.   In  this  case,  the  composite 
estimate  is: 

x  =  (Zxi/ai)/(El/ai)  ,  y  =  {ly±/o±) / (ll/o±)     , 

2  2      2  2 

a-   =  n/(Zl/a.)   ,  a~   =  n/(El/a.)   and  o~~    =    0. 
x  i      y  l        xy 

In  this  example,  since  x  and  y  are  determined  by  a  circular 
normal  distribution,  the  minimum  area  confidence  regions  are 
circles  and  orientation  is  not  an  issue. 
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y  North 


Figure  5.   The  ellipses  represent  the  confidence  ellipses  of 
the  first  example.   The  position  estimates  are  at  the  center 
of  the  ellipses.   The  numbers  indicate  the  order  of  the  esti- 
mate in  the  table  on  Page  17. 


19 


References 

1.  Bowditch,  N. ,  American  Practical  Navigator,  Vol.  I,  Appendix 
Q,  Defense  Mapping  Agency  Hydrographic  Center,  1977. 

2.  Daniels,  H.  E.,  "The  Theory  of  Position  Finding,"  J.  Royal 
Stat.  Soc.  (B) ,  Vol  13,  1951,  pp.  106-207. 

3.  Forrest,  R.  N.,  "Programs  for  a  Target  Position  Estimation 
Procedure,"  NPS71-83-002 ,  Naval  Postgraduate  School,  Monterey, 
CA  93943,  March  1933.* 

*Note:  The  angle  labeled  y  in  Figure  4  of  this  report  is 
the  ellipse  major  axis  direction.  The  angle  between  the 
y  and  y1  axis  should  be  labeled  y  for  the  case  indicated 
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